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Precise measurement of the Higgs boson properties is an important issue of the Interna- 
tional Linear Collider (ILC) experiment. We studied the accuracy of the Higgs mass re- 
construction in the ZH — > qqH multi-jet process with the Higgs mass of Mh — 120 GeV 
at y/s — 250 GeV with the ILD detector model. In this study, we obtained the recon- 
structed Higgs mass of Mh = 120.79 ± 0.089 GeV and 5.3% measurement accuracy of 
the cross-section for ZH — > qqbb with the integrated luminosity of C — 250 fb" 1 data 
samples. 



1 Introduction 



International Linear Collider (ILC) pQ is a future e + e~ collider experiment for the precise 
measurement and the validation of the Standard Model (SM) physics, especially for the 
measurement of the Higgs boson property, even the discovery of the Higgs boson will be 
realized in Large Hadron Collider (LHC) experiment. In the SM, light Higgs boson mass 
(M H ) is predicted around the 114.4 GeV < M H < 160 GeV from the study in LEP g] 
and Tevatron [3] experiment. The largest production cross-section for SM Higgs boson is 

-> ZH) process which associated with 
as shown in Fig. [IJ around the ZH 



obtained through the Higgs-strahlung (e + e~ 
the Z boson and the Z mainly decays to qq pair, 
production threshold energy shown in Fig. [2] (a). 

Since Higgs boson mainly decays to bb pair at 
the Higgs mass below 140 GeV region as shown in 
Fig. [2] (b), the final state of the ZH — > qqH process 
forms the four-jet. In ILC experiment, the most of 
interesting physics processes including ZH process 
form the multi-jets final state from the decay of gage 
bosons (W,Z) and heavy flavor quarks (b,c), thus 
ILC detectors are required to have the good jet en- 
ergy resolution for the precise measurement. There 
are three detector concepts, SiD, ILD and 4 th for the 
ILC detector, and ILD is the merged concepts of the 
previous GLD |U (Asian group) and LDC [S] (Eu- 
ropean group) models for the Letter of Intent (LOI) 
submission [6|. In order to achieve the best jet en- 
ergy resolution, ILD adopt the Particle Flow Algo- 
rithm (PFA) suited detector design. Since the PFA 

performance is degraded by the cluster overlapping and the double-counting of the particles 
energy in the calorimeter, particles separation in the calorimeter is an important key for 
better PFA performance. 




Figure 1 : Higgs boson production via 
Higgs-strahlung (ZH) process and Z 
mainly decay to qq. 
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Figure 2: (a). Production cross-section of the Higgs boson as a function of center-of-mass 
energy (y/s) and (b). branching ratio of the Higgs decay as a function of the Higgs mass. 



The figure-of-merit of the PFA performance from each detector parameter relating to the 
particles separation in the calorimeter is described as F.O.M. = BR 2 / y/ a 2 + Rm 2 , where B 
is a magnetic field, R is a detector radius, a is a segmentation of the calorimeter and Rm is 
a effective Moliere radius of the calorimeter. In order to maximize the F.O.M. , ILD detector 
adopts the large radius tracker and high granularity calorimeter with 3.5 T magnetic field. 
In this analysis, we study the direct reconstruction of the Higgs boson mass with the full 
detector simulation for ZH — > qqH , H — > bb four-jet mode with the ILD detector model. 

2 Simulation tools 

For full detector simulation study, we use the ILD detector model based Monte Carlo 
(MC) full simulation package called Mokka, which is based on the MC simulation pack- 
age Geant4 [7]. Generated MC hits are reconstructed and smeared in the reconstruction 
package called MarineReco which includes the PFA package called PandoraPFA [8]. Since 
y/s = 250 GeV reconstructed and skimmed signal and background samples called DST files 
are generated for the LOI physics analysis in ILD group, we use these DST data samples 
saved in the linear collider common data format called LCIO. For the DST data sample 
analysis, we use the useful analysis package library called Anlib for the event shape analysis 
and jets reconstruction, and analysis process is handled through the Root [9] based analysis 
framework called JSF [TU]. For the comparison of the PFA performance between realistic 
PFA and perfect-clustering PFA, we also use the GLD detector model MC full simulator 
called Jupiter [TT] with the generating the signal and background events by PYTHIA, and 
reconstruction package called Satellites [12] based on Root, both of them are also controlled 
in the JSF framework. From the comparison of the ZH — > qqH in GLD detector model, 
shown in Fig [21 PandoraPFA reconstruction performance (a) achieve the comparable perfor- 
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mance with perfect-clustering PFA (b) in terms of the reconstructed Higgs mass distribution 
width of a which corresponds to the jet energy resolution even only the ZZ — > qqq'q' back- 
ground is considered. Therefore, we shift to the full SM background analysis with common 
DST data. 




Higgs mass (GeV) Higgs mass (GeV) 

(a) (b) 

Figure 3: Comparison of the reconstructed Higgs mass distribution for ZH — > qqbb only with 
ZZ background in GLD detector model with the different PFA clustering of (a) realistic 
PandoraPFA and (b) perfect clustering PFA. 



3 Analysis Procedure of ZH — > qqH mode 
3.1 MC samples 



The SM Higgs boson is mainly produced through 
the Higgs-strahlung e + e~ — > ZH process around the 
production threshold center-of-mass energy {^fs ~ 
230 GeV). Since the main decay mode at Mjj < 2,Mw, 
Higgs boson mainly decays to bb pair, thus largest pro- 
duction cross-section is obtained from the ZH — > qqbb 
process, which forms four-jet final state and both Z 
and H can be reconstructed directly. Fig. U] shows 
the typical event display of the ZH — ► qqH in JSF. 
In this analysis, we assume the center-of-mass energy 
as the ZH production threshold of y/s = 250 GeV 
and the light Higgs mass of Mh = 120 GeV. Each 
DST data samples is scaled to the integrated lu- 
minosity of C = 250 fb _1 and the beam polariza- 
tion to P(e+,e _ ) = (30%, -80%). The main back- 
grounds for ZH — > qqbb are considered as following 
processes: ZH -> -> qq, e+e~ -> WW/ZZ -> 




Figure 4: Typical event display of 
the ZH — > qqH four-jet final state. 
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qq'q"q"' or qqq'q', e + e — > WW — > viqq' and e + e — * ZZ — > Generated signal and 

background MC samples which scaled to be C = 250 fb -1 are summarized in Table. [T] 



MC samples (C = 250 ftr 1 ) 


ZH -> OT F (sig) 


qqqq 


viqq 


me 


qq 


Number of generated events 


51763 


814163 


302807 


98127 


2529928 



Table 1 : Generated signal and background MC data samples scaled with L = 250 fb 1 . 

In order to correct the escape energy from the heavy quark decay including neutrinos, 
kinematic five constraint (5C) fit is applied, which consists of the four constraints (4C) of 
momentum balance P x ,y,ZA = 0) and jets energy balance (^2 Ej — y/s = 0) of the four-jet 
and one Z mass constraints for Z candidate di-jet. For the kinematic fitting, jet energies 
(Ej) and jet angles (6, 0) of each jet are used as measured variables. Finally, reconstructed 
Higgs mass distribution is fitted with the Gaussian convoluted with Gaussian function for the 
signal and exponential function for the contribution from background events which remain 
after the Higgs boson selections. 

3.2 Jet Reconstruction 

Since the final state of the ZH — > qqH mode forms four-jet, after the PandoraPFA clustering, 
forced four-jet clustering based on Durham jet-clustering algorithm has applied. In order to 
select the best jet pair combination from the four-jet, following x 2 value is evaluated, 

2 = ( My - Mz V { Misshhi-M z 

V a M z ) V &MM H 

where M12 is Z candidate di-jet mass, MissM^ is a missing mass of the remaining Higgs 
candidate di-jet, Mz is the Z boson mass (91.2 GeV), and om z and <JMi SS M 3l are sigma of 
distribution of the reconstructed Z boson mass and the missing mass of the Higgs candidate 
jets, respectively. In order to select the best jets pair combination, x 2 < 10 is required for 
the reconstructed jets pair. 

3.3 Event selection 

After the x 2 cut to select the best jet pair combination, following event selections are applied 
for background rejection: 

(a) visible energy : 200 < E ms < 270 GeV; 

(b) Longitudinal momentum of the Z : \Pg z \ <70 GeV to reduce ZZ background; 

(c) Higgs production angle : | cos 9h\ < 0.85 to reduce the ZZ background; 

(d) thrust angle : thrust < 0.9; 

(e) Number of particles: N part i C ie > 40 to suppress the till background; 

(f) Maximum and minimum jet energy fraction: E m i n j E max > 0.25; 

(g) Maximum momentum of jet: Pa < 100 GeV; 

\°/ J J max ' 

(h) Y Plus : YPlus > 0.0001; 
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(i) Y Minus : YMinus > 0.001; 

(j) Minimum angle of Z-H jets : 20 < 9 Z Hj min < 135; 
(k) Maximum angle of Z-H jets : 110 < 0zHj max ', 
(1) 6-tagging : P btag > 0.5 from LCFIVTX package. 

The distribution and its cut positions for each selection variable are shown in Fig. [5j 
Since the W/Z generated in the WW/ZZ background event are relatively boosted compare 
to the Z generated in ZH signal event, longitudinal momentum of Z (Pez) an d maximum 
momentum in jets {Pj max ) are higher in WW/ZZ background event than in signal event. 
None jet-like background events are reduced by the number of particles (Nppo) cut. Y 
Plus and Y Minus values are threshold Y-values used in the jet clustering topology which 
reconstructed from four-jet to five-jet or three-jet, respectively. Minimum and maximum 
angles between Z and H candidate jets are also used for the separation by the event shape 
difference between ZH event and backgrounds. 



■70<zpl&&zpl<70 Maqh 200<evisSSevis<270 Mnqh pjmax<100 Mqnh hj1yplus>1.0e-4fl&hj2yplus>1.0e-4 




# of panicles E Jet max/min ratio z-H Jets mln angle B tagging probability 



Figure 5: Distribution of each selection variable and its cut positions to select ZH — ► qqbb 
event. 

Finally, we apply the vertex tagging selection for the neural net output of the 6-likeness 
analyzed in the vertexing package called LCFIVTX in ilesoft. The reduction summary in 
each event selection is listed in the Table [21 
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Selections 


ZH -> qqH(Sig) 


9999 




UU 


99 


no cuts 


51745 


814162 


302807 


98127 


2529928 


X 


^fi74H (7~\ D9 °7r,\ 

JUI^O ^1 /o j 


688703 (84 5Q %'! 


1 Q043 (F> 9Q %"! 


95375 (95 86 %\ 


541 859 (91 49 %'! 


\ jr iz\ 




479403 f 58 88 %1 

T: 1 i7TUO I UU.OC ' ^ / 


1 9832 (A 24 %1 


5565 (5 67 %1 


993883 M 1 69 %1 


E vis 


34924 (67.49 %) 


477994 (58.71 %) 


12457 (4.11 %) 


5335 (5.44 %) 


287324 (11.36 %) 


| COS &H 


30451 (58.85 %) 


397270 (48.79 %) 


9934 (3.28 %) 


2167 (2.21 %) 


223873 (8.85 %) 


thrust 


29916 (57.81 %) 


389703 (47.87 %) 


8312 (2.75 %) 


1422 (1.45 %) 


103283 (4.08 %) 


NparticUs 


29820 (57.63 %) 


389514 (47.84 %) 


4353 (1.44 %) 


(0.00 %) 


87022 (3.44 %) 


Ei . /Ej 


27843 (53.81 %) 


297580 (36.55 %) 


1603 (0.53 %) 


(0.00 %) 


40880 (1.62 %) 


mn.r. 


27622 (53.38 %) 


289490 (35.56 %) 


1500 (0.50 %) 


(0.00 %) 


31382 (1.24 %) 


Yplus 


27607 (53.35 %) 


288421 (35.43 %) 


1465 (0.48 %) 


(0.00 %) 


30773 (1.22 %) 


Yminus 


27559 (53.26 %) 


287825 (35.35 %) 


1354 (0.45 %) 


(0.00 %) 


27250 (1.08 %) 


^Z — Hj „,j n 


27311 (52.78 %) 


285704 (35.09 %) 


1284 (0.42 %) 


(0.00 %) 


24601 (0.97 %) 




27031 (52.24 %) 


277203 (34.05 %) 


1263 (0.42 %) 


(0.00 %) 


24280 (0.96 %) 


b — tagging 


5972 (11.54 %) 


4732 (0.58 %) 


(0.00 %) 


(0.00 %) 


458 (0.02 %) 



Tabic 2: Backgrounds reduction summary in each selection for ZH — > qq&fr. 

From the reduction summary of Table. ££££ four-leptonic background can be suppressed 
completely by number of particles cut (Nppo s < 40) and the remaining backgrounds are 
qqqq and qq which including 6-quarks event after applying the 6-tagging. 

4 Results 

Reconstructed Higgs mass distribution after the se- 
lection of ZH — > qqbb is fitted with the function 
of Gaussian convoluted Gaussian with the exponen- 
tial function assuming the background, as shown in 
Fig. [51 Fitted results of the reconstructed ZH — > 
qqbb Higgs mass distribution are summarized in the 
Table. [31 From the fitted results, Higgs mass (Mh = 
120 GeV at MC) is reconstructed as M H = 120.79 ± 
0.089GeV and the measurement accuracy of cross- 
section to ZH — > qqbb is obtained as 5a ja = 5.3%. 

5 Conclusion 



Simulation study of the direct reconstruction of the 
Higgs boson in ZH — ► qqbb four-jet mode with the 
Higgs mass of 120 GeV at the ^/s = 250 GeV and Figure 6: Reconstructed Higgs mass 
the integrated luminosity of C = 250/6 -1 has per- distribution of ZH — > qqbb. 
formed for the ILD detector model considering with 



Higgs mass (M H = 120 GeV at MC) 


M H = 120.79 (GeV) 


Measurement accuracy of Mh 


SM H = 89 (McV) 


Measurement accuracy of a(ZH — > qqbb) 


Sa/a = 5.3% 



Table 3: Fitted results for the reconstructed Higgs mass distribution. 
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the qqqq, vtqq ££££, qq background processes. From 

the study, measurement accuracy of the reconstructed Higgs mass is estimated as 87 MeV 
and the measurement accuracy of the cross-section of ZH — > qqbb mode is obtained as 
Safer = 5.3%. 
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